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SUMMARY 

The solution of heat-transfer problems has become vital for many 
aeronautical applications. The shapes of objects to be cooled can 
often be approximated by cylinders of various cross sections with flow 
normal to the axis as, for instance, heat transfer on gas-turbine blades 
and on airfoils heated for deicing purposes. A laminar region always 
exists near the stagnation point of such objects. 

A method previously presented by E. R. G- Eckert permits the calcu- 
lation of local heat transfer around the periphery of cylinders of arbi- 
trary cross section in the laminar region for flow of a fluid with con- 
stant property values with an accuracy sufficient for engineering 
purposes. The method is based on exact solutions of the boundary-layer 
equations for incompressible wedge-type flow and on the postulate that 
on any location of the cylinder the boundary-layer growth is the same as 
that on a wedge with comparable flow conditions. This method is extended 
herein to take into account the influence of large temperature differ- 
ences between the cylinder wall and the flow as well as the influence of 
transpiration cooling when the same medium as in the outside flow is used 
as coolant. Prepared charts make the calculation procedure very rapid. 
For cylinders with solid walls and elliptic cross section, a comparison 
is made between the results of calculations based on the presented method 
and the results of calculations by other known methods as well as of 
experimental investigations. 


INTRODUCTION 

A calculation of heat transfer to cylinders with arbitrary cross 
section in an air flow normal to the axis by a solution of the boundary- 
layer equations is a difficult problem, even for the laminar region. The 
problem is especially complicated by the large number of parameters 


2 


NACA TN 2733 


influencing heat transfer. Such parameters are the shape of the cross 
section of the cylinder, the Mach number which determines the flow out- 
side the boundary layer, the temperatures on the surface of the cylinder 
as well as in the stream, the stream velocity determining the internal 
heat generation, and the temperature distribution around the circumfer- 
ence of the cylinder . If the cylinder is cooled by the transpiration- 
cooling method in which a coolant is ejected through the porous surface 
into the outside stream, the amount of coolant and its distribution 
around the circumference of the cross section of the cylinder are addi- 
tional parameters. Even if a solution is obtained for such a problem, 
for instance by use of an electronic computer, this solution is very 
restricted because of the many parameters. Up to the present time, 
therefore, the problem has been attacked only under simplifying 
restrictions. 

The restrictions most commonly used are: (l) low velocities, (2) 

constant property values, (3) constant wall temperatures, and (4) imper- 
meable surfaces (no transpiration cooling). Under restriction (2), the 
development of the boundary layer along the cylindrical surface is inde- 
pendent of the heat transfer} available knowledge on the flow boundary 
layer can therefore be used as a basis for a heat -transfer calculation. 
Under the simplifying assumptions, which are necessary in order to 
transform the general viscous -flow equations into the boundary-layer 
equations, the development of the flow boundary layer does not depend 
immediately on the shape of the cross section of the cylinder, but only 
on the velocity distribution in the stream outside the boundary layer 
and along its surface. 

One method which was applied successfully to obtain a solution of 
the flow boundary- layer equation developed the stream velocity along the 
surface of the cylinder in a power series of the distance from the stag- 
nation point measured along the circumference of the cylinder. In ref- 
erence 1 this method is used to solve the heat— transfer problem. It is 
also shown that the temperature field within the boundary layer can also 
be presented in a power series of the distance from the stagnation point 
in which the single terms contain only universal functions of a dimen- 
sionless wall distance and of the Prandtl number of the fluid. The heat 
transfer to the surface is given by an analogous series with terms 
depending on the Prandtl number* The calculation of the universal func- 
tions, however, is a tedious process, and accordingly these functions 
are known only for a limited number of terms. For air with a Prandtl 
number of 0.7, they are presented in reference 1. For a gas with a 
Prandtl number of 1, they are contained in reference 2, which is based 
on reference 3, in which the boundary -layer flow of a yawed cylinder is 
calculated. The fact that the boundary -layer equation for the velocity 
component parallel to the axis of a yawed cylinder is identical in form 
to the boundary-layer equation describing the temperature field for a 
fluid with a Prandtl number of 1, flowing normal to the axis of the 
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cylinder, was used in reference 2 to determine heat transfer to such 
cylinders. The presentation of more terms of the series is announced in 
reference 4. It is found, however, that the velocity distribution for 
only a limited range of cross sections of cylinders can he represented 
by a power series converging rapidly enough that the number of the known 
universal functions is sufficient to calculate the heat transfer. 

The difficulties connected with a solution of the boundary-layer 
equations point out the need for an approximate approach which, with a 
small expenditure of time, would determine heat-transfer coefficients 
with an accuracy sufficient for engineering purposes. A considerable 
number of such approaches were tried in the past with results which 
differ greatly as shown in figure 1, taken from reference 2. 

The simplest procedure is probably the use of the heat-transfer 
coefficients as calculated in reference 5 for flow along a flat plate 
with a constant velocity. The fact that in reality the stream velocity 
varies along the cross section of the cylinder is taken into account by 
calculating the local heat-transfer coefficients with the velocity found 
in the stream at the considered distance from the stagnation point. 

This method is contained in a summary presented in reference 6. Unfor- 
tunately, such an approach gives heat-transfer coefficients which are 
considerably low in many cases (see fig. 1). 

Better agreement was obtained by another approach (reference 7) 
which uses, instead of the flat -plate solution, a family of solutions of 
the boundary -layer equations which can be obtained in a general form, 
namely, for the case where the stream velocity varies along the surface 
as a certain power of the distance from the stagnation point. Such a 
velocity vaxiation is obtained in incompressible flow around wedges. 

The solutions for such a type of flow were used to obtain approximate 
heat-transfer coefficients for a cylinder with arbitrary cross section 
by stipulating that the local heat-transfer coefficient on any location 
along the cylinder is identical to the local heat-transfer coefficient 
on a wedge for which, at the same distance from the stagnation point, 
the stream velocity and its gradient are the same as those on the inves- 
tigated cylinder. This approach was subsequently used by different 
authors, and is described, for instance, in references 8 and 9. It 
takes into account the stream conditions which influence the boundary- 
layer growth at the location at which the heat transfer is going to be 
determined; however, it does not properly account for the development of 
the boundary layer in the range upstream of the point considered. This 
development may be different on the cylinder and on the equivalent wedge. 

Another group uses an integrated momentum equation for the boundary- 
layer flow as proposed by von Karman and K- Pohlhausen (references 10 
and 11) to calculate the velocity boundary layer. Different procedures 
were proposed to determine local heat-transfer coefficients from the 
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known velocity boundary layer. Some investigators use Reynolds analogy 
directly (reference 12) or with a correction for Prandtl numbers diiferent 
from 1 (reference 13). Such approaches give heat-transfer coefficients 
which are considerably high in many cases , as shown, for instance, in fig- 
ure 1. More accurate results were obtained when the heat transfer was 
determined by solving an integrated heat-flow equation for the boundary 
layer. The velocity field within the boundary layer has to be known in 
this approach, since the flow velocities within the boundary layer occur 
in the mentioned heat -flow equation. This method was originated by 
Kroujiline (reference 14) . Extensions and simplifications are contained 
in references 15 to 18 and an extension to compressible flow of a fluid 
with a Prandtl number equal to 1 is found in references 19 and 20. Use- 
ful information is also contained in a summarizing report (reference 21) . 

Another approach starts from a consideration of the fact that the 
use of the heat-transfer coefficients for wedge-type profiles as 
described previously was found to give heat-transfer coefficients with 
a fairly good accuracy. It should be expected that these heat-transfer 
coefficients can be improved to a degree which is sufficient for all 
engineering purposes by a method which takes into account in some approx- 
imate way the previous history of the boundary layer- Such a method, 
called the equivalent wedge-type flow method, was proposed in refer- 
ence 22, extended to heat transfer at high flow velocities and variable 
wall temperature in reference 23, and extended to transpiration cooling 
with small temperature differences in reference 24. The advantages of 
this method are that no knowledge of the velocity boundary layer is 
required and that it can be readily extended to take into account the 
effects of large temperature differences, of transpiration cooling, and 
of variable wall temperature as soon as the corresponding solutions for 
the wedge -type flow are available. 

This report presents such an extension, made at the NACA Lewis lab- 
oratory, which includes the effects of large temperature differences and 
of transpiration cooling. It is based on exact boundary -layer solutions 
for wedge -type flow with large temperature differences and with trans- 
spiration cooling (reference 25) . Charts were prepared which make the 
calculation of heat transfer around cylinders of any arbitrary cross 
section more rapid. 


SOLUTION OF BOUNDARY -LAYER EQUATIONS FOR WEDGE -TYPE FLOW 
Boundary-Layer Equations 

The following boundary-layer equations describe the velocity and 
temperature fields in a laminar steady two-dimensional gas flow: the 

momentum equation, the continuity equation, and the energy equation. 
The momentum equation is 
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, Bu Bu B / Bu\ Bp 
pu BI + pv B^ “ B^F St) ~ 5£ ^ 

when body forces are neglected. (All symbols are defi n ed in appendix A) 
consistent units are used throughout the report.) Since the pressure 
variation normal to the surface throughout the boundary layer may be 
neglected, it follows that the pressure is prescribed by the conditions 
in the stream outside the boundary layer and can be connected with the 
velocity u s in the stream and just outside the boundary layer by the 
Bernoulli equation 


Bp du s 

~ ** p s u s 3sr 

The introduction of this expression changes the momentum equation to the 
form 

Bu Bu B / Bu) ^ u s , s 

p Bx + pv By “ By Y* By/ + p s U s Sx~ (2) 

The continuity equation is 

^ (pu) + ^ (pv) » 0 (3) 

and the energy equation is 



The heat generated by internal friction, described by the second term on 
the right side of equation (4), and the temperature variation connected 
with expansion, described by the third term, can be neglected as long as 
the difference between the total and the static temperature in the gas 
stream is small compared with the difference between the wall temperature 
and the temperature in the gas stream. For this condition, then, only 
the first term on the right side of equation (4) is retained, and the 
energy equation assumes the form 

/ Bt Bt\ B /, Bt) 

pC P \ U Bx + v By/ By \ By/ 

Equations (2), (3), and (5) include the case of transpiration cool- 
ing when the same medium as that in the outside flow is used as coolant 
and the boundary conditions are properly defined. 
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u s 0, v at v w , and T = T w when y = 0 
u -) u s and T T s when y 00 


( 6 ) 


The property values p., k, Cp, and p appearing in the equations depend 
on temperature and pressure. The variation with pressure can he neg- 
lected at the small velocities to which the energy equation was already 
restricted by disregarding the internal friction and the expansion terms. 
The influence of the temperature dependency, however, may be appreciable 
in applications with large temperature differences within the boundary 
layer. Solutions of the boundary -layer equations which take into account 
the temperature variation of the property values were obtained in ref- 
erences 9 and 25, in which the partial differential equations were 
transformed into total differential equations. 


Change of Variables 

The transformation of the partial differential equations into total 
differential equations is possible under the following specialized con- 
ditions: The stream velocity is assumed to vary as a power function of 

the distance from the stagnation point measured along the surface of the 
cylinder . 


u s H (7) 

It has recently become customary to refer to the exponent m in this 
equation as "Euler number." The Euler number can be expressed by the 
Bernoulli equation in the following way: 


m 


dp 

' 33E 


P u 
H s s 


X 


( 8 ) 


In addition, the temperature of the wall is assumed to be constant and 
the property values are assumed to vary proportionally to a power of 
the absolute temperature T. The numerical calculations were made for 
air. The exponents used were 0.7 for the viscosity, 0.85 for the heat 
conductivity, 0.19 for the specific heat, and -1*0 for the density. 

The variables 
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(9) 
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are used to transform equations (2), (3), and (5) into total differen- 
tial equations presenting f and Q as functions of rj only. The 
stream function f appearing in equations (9) is defined in such a way 
as to eliminate the continuity equation (3) . 


pu 


pv 


d(o w t) 

dy 

d(p w *) 

5x 




( 10 ) 


Introducing the new variables into the second of equations (10) leads 
to the following expression for the velocity component normal to the 
surface : 





The velocity at the surface itself follows: 


(ID 


v w " ~ 


1+m 


w 


' ^w u s 

Pw* 


( 12 ) 


The transformation therefore prescribes a certain variation of the cool- 
ant velocity v w along the surface, since the function f w has to be 
constant (independent of x) . Since the stream velocity is described 
by equation (7), the coolant velocity v w is also proportional to some 
power of x. it is shown in reference 26 that such a variation of the 
coolant velocity leads to a constant wall temperature and is therefore 
consistent with the assumed constant wall temperature when heat transfer 
by radiation may be neglected. The transformed equations are presented 
in references 9 and 25, together with the solutions for a Prandtl num- 
ber Pr of 0.7, and for a range of Euler number m, temperature ratio 
T s /t w , an( l Ibe parameter f w describing the cooling-air flow through 
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the porous surface. The results contain expressions for the thickness 
of the flow boundary layer which are defined in two ways: the dis- 

placement thickness 



jpu_\ 

Ps U S,/ 


<iy 


and the momentum thickness 


&i = 


Pop 


Jo 


pH 

Ps u s 



(13) 


(14) 


The thermal boundary layer is characterized in this report by the con- 
vection thickness 


'00 


6c - 


T - T 
pu s 


p s u s ^w " 


ay 


(15) 


In addition, a thermal boundary-layer thickness will be used herein and 
is defined in the following way: 


8 t ■ 




(16) 


Values for this boundary-layer thickness can be easily calculated from 
results presented in reference 25. 


Application to High Velocities 

The solutions described apply exactly only to flow with small 
velocities. Practically, the limiting velocity up to which it is pos- 
sible to neglect the frictional and the expansion terms can be set quite 
high for a gasj this fact can be understood from the following trans- 
formation of the energy equation, in which only the specific heat is 
regarded constant. If the momentum equation (l) is multiplied by the 
velocity u and added to the energy equation (4) and if, in addition, 
the total temperature Tqi = T + u^/2cp is introduced, the following 

expression is obtained: 


P c p 


3^ 


St, 


+ v 


T 



( k *t) 

l s 

76 

i — i 
r 

L 

/By 

\ k 3T/ 

+ 2 37 

\S V- ‘ Er/ 


(17) 
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The last term on the right side of the equation vanishes for a Prandtl 
number equal to 1. In this case, the energy equation has the same form 
as the one for low velocities in which the friction and the expansion 
terms were neglected. The only difference lies in the fact that the 
total temperature appears in the energy equation. For cases which have 
a Prandtl number not too far from 1, the last term in equation (17) will 
be comparatively small up to considerable velocities, and the energy 
equation (5) used in the following considerations applies to this condi- 
tion when the temperature T is interpreted as total temperature. It 
will be shown later that as far as heat transfer is concerned, the range 
in which the results of a calculation with equation (5) may be used can 
be extended even further by using a properly defined adiabatic wall tem- 
perature instead of the total gas temperature. 

The property values p, k, Cp, and p depend for gases on the tem- 
perature. This dependency was taken into account in the described cal- 
culations. The density depends, in addition, on the pressure, and the 
pressure variation may become considerable at high Mach numbers. There 
are indications, however, that calculations which neglect this pressure 
variation can be used with sufficient accuracy over the whole subsonic 
range, as is pointed out in reference 27 by an investigation of results 
obtained by L. Howarth (reference 28). 


EXTENSION OF THEORY TO ARBITRARY BODIES 
Determination of Equivalent Wedge 

The solutions discussed in the previous paragraph are in an exact 
sense restricted to a certain type of velocity variation along the 
cylindrical surface, namely, a stream velocity which just outside the 
boundary layer is proportional to some power of the distance from the 
stagnation point. Such a velocity distribution is realized, for 
instance, in incompressible flow around wedges. The wedge-type solu- 
tions may be used, however, to obtain approximate heat-transfer coeffi- 
cients on cylinders of arbitrary cross section. One approach in this 
direction assumes that the heat-transfer coefficient on any point along 
the circumference of a profile with arbitrary cross section is the same 
as that on a wedge at the same distance from the leading edge , provided 
the stream velocity and its gradient on the wedge and on the arbitrary 
profile have the same value at the location considered and that the tem- 
perature ratio T s /T w is the same . It will be shown that such an 

approach takes into account the right stream conditions at the local 
spot for which the heat-transfer coefficient is to be determined. How- 
ever, it does not properly consider the previous history within the 
boundary layer. Numerical calculations presented later show that heat- 
transfer coefficients obtained in such a way are in most cases within 


10 


NACA TN 2733 


about 15-percent agreement with experimental data. It is to be expected 
that a modification which accounts in some approximate manner for the 
conditions in the boundary layer upstream of the point under considera- 
tion should improve this approximation to the desired degree. This mod- 
ification is made in reference 22 by the stipulation that the rate of 
increase of the boundary-layer thickness is the same on the considered 
point of an arbitrary profile and on the point of a wedge which has 
the same boundary-layer thickness, the same stream velocity, and the 
same stream velocity gradient. This same stipulation will be used in 
the present report. For a given temperature ratio Tg/T-w, the heat- 
transfer coefficients on a wedge depend on the Euler number m and the 
value f w characterizing the coolant flow through the porous surface. 
These parameters which define the equivalent wedge profile will now be 
expressed by the boundary-layer thickness and the local stream velocity 
gradient . 

For the wedge-type profile, the stream velocity is expressed by 
the power law 

u s = Cf 1 (18) 

in which the value i expresses the distance from the leading edge 
measured along the wedge surface in order to distinguish it from the 
distance of the point under consideration from the stagnation point on 
the arbitrary profile, which will be denoted by x.' The variables used 
for the transformation of the original boundary-layer equations in the 
previous section may now be written 


T] = 


Pw u s 




and 


v A EL 

m+1 w yPw u s 


(19) 


( 20 ) 


r 


To a certain value y indicating the boundary-layer thickness 6 
belongs a value % of the coordinate T) defined by the equation 


% 


Pw u s 


( 21 ) 
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In order to eliminate the distance £ from this equation, equation (18) 
is differentiated to obtain 

Su c .. u c 

= xnC| m_1 = m -|- ( 22 ) 

Since the velocity gradient on the wedge profile is assumed the same as 
that on the profile under consideration, it follows that Su s /S £ = du s /chc. 
This equality gives for the coordinate | the expression 


du s / dx 


(25) 


When this egression is introduced into equation (2l), there is obtained 

_ N / du 

\dx 


In this expression, % (denoted as (6/x) '^Re in reference 25) is a 
function of the Euler number m and of the coolant-flow parameter f W* 

If this equation is therefore written in the form 

2 P w 6 2 dUg 

^b 111 |a w dx 

the left side is a function of m and f v and equation (24) relates 
both values to the boundary-layer thickness 5 and the velocity gradi- 
ent du s /dx. In order to obtain a second relation for m and f v , the 
coordinate £ is replaced in equation (20). The result is 

m+l ^ Pv v v & 

- -r f »% - — 

which is written again in such a way that the left side is a function of 
the Euler number m and the flow parameter fw> which can be calculated 
from the results in reference 25. Both equations (24) and (25) are 
therefore sufficient to determine the equivalent wedge profile. 
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Equations for Boundary-Layer Thickness and Heat Transfer 


The next step is to develop a differential equation for the 
houndary-layer thickness from the postulate that the houndary-layer 
gradient d8/dx is the same for the real profile as for the equivalent 
wedge profile. For the wedge profile, the houndary-layer thickness is 
given hy the expression 



which is obtained from equation (2l) hy replacing the stream velocity 
with equation (18) and solving for the houndary-layer thickness. A dif- 
ferentiation of this equation and the use of equations (23) and (24) 
result in 


d5 _ 1-m du s/^ _ 1-m 2 ^ 

dx 2m u s 2 h p^-SUg 


(26) 


This is a differential equation for the houndary-layer thickness which 
contains only values which are lmown for the profile under consideration 
or which are determined from equations (24) and (25) for the equivalent 
wedge-type flow. An integration of the differential equation gives the 
boundary layer along the circumference of the profile under consideration. 

The local heat -transfer coefficient is defined hy the following 
equation: 


h(T s -T w ) 



Introducing the dimensionless temperature ratio given in equation (9) 
and the coordinate £ results in 

* - w I - w ^ < 27 > 

The heat-transfer coefficient may he calculated from this expression as 
soon as the houndary-layer thickness 6 is known, since and T] b 

are functions of m and f w contained in reference 25. 
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Up to the present time, no recommendation has been made as to 
which boundary -layer thickness should be used in the prescribed proce- 
dure. When the momentum thickness is used in the foregoing equations, 
it is easily understandable that the integrated momentum equation is 
satisfied and the method of calculation becomes the same as the one pro- 
posed by von Karman in reference 10. This fact can be proved mathemati- 
cally by a procedure completely analogous to the one used in appendix B. 
On the other hand, the use of the convection thickness as defined in 
equation (15) satisfies the integrated heat-flow equation within the 
boundary layer, as shown in appendix B. The use of both boundary- layer 
thicknesses leads to somewhat different results for the local heat- 
transfer coefficient, and the question arises as to which is preferable. 
It is pointed out by Schuh in reference 23 that for the purpose of deter- 
mining heat-transfer coefficients it is more important to satisfy the 
heat-flow balance, and the use of the convection thickness was therefore 
recommended. In reference 22, the use of the thermal boundary-layer 
thickness as defined in equation (16) is investigated, and the results 
of the calculation with this boundary-layer thickness are found to 
agree even better with measured values and with other calculations. The 
convection thickness 6 C and the thermal thickness 6^ for the boun- 
dary layer will therefore be used in parallel in the following numerical 
evaluations. 


CALCULATION PROCEDURE 


Use of Dimensionless Variables 

The procedure which may be followed in determining local heat- 
transfer coefficients with the relations developed in the preceding sec- 
tion is now explained. Figure 2 shows a sketch of a cylinder with arbi- 
trary cross section and the notation used in the analysis. Before 
numerical calculations are made, however, it is advisable to change to 
dimensionless quantities. In order to make this change, the distance x 
is divided by the major axis L of the cylinder and the mass velocity 
in the direction of x is divided by an upstream mass velocity. All 
lengths and mass velocities parallel to y are, in addition, multiplied 
by the square root of the Reynolds number Reg based on the major axis 
and the upstream mass velocity. The dimensionless variables which are 
subsequently needed are 



(28) 


5 * = 



(29) 
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where 


* Pw u s 

u =" SJT 


"0 s,0 


P v 


Vw Pn u , 


w w 




0 s,0 


0 


Be 0 = ^f0 
0 


(30) 


(31) 


(32) 


By use of these dimensionless quantities, equation (26) is transformed 
into 


where 


d6* M 
dx + “ S*u* 


M = 


1-m 

2 


% 


= M 


du 

s 

dp 


8*1 


★ 

v w 6 


■) 


(33) 


(34) 


according to equations (24) and (25), which, in dimensionless values, 
are 




and 


m+1 

~ 


•^w% - 


Introduction of the dimensionless quantities into equation (27) 
leads to 


Nu = H_ 

&* 


where 


and 



n = 




(35) 


(36) 


(37) 
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Charts and Calculation Procedure for Prescribed Coolant Flow 

Charts have been prepared which present the functions M and N as 
expressed by equations (34) and (37) in dependence on (dus/dx*)5* 2 and 
v*S*. When the values presented in reference 25 were considered for the 
construction of these charts, it was found that in certain cases one of 
the free-stream boundary conditions was not fulfilled. This same condi- 
tion had been overlooked by other previous investigators . As a conse- 
quence, spot r ecal culations were made; it was found that the maximum 
error in Nu/ V^O resulting from these recalculations was of the order 
of 6 percent. Consequently, no attempt was made to recalculate the entire 
field because of the large amount of labor involved. The charts presented 
herein were constructed from a combination of these recalculations and 
results presented in reference 25. Even though all recalculations were 
not made, it is believed that the charts as presented will give results 
with an accuracy of the order of 2 or 3 percent. In figures 3 and 4, the 
dimensionless convection thickness of the boundary layer is used; in fig- 
ures 5 and 6, the dimensionless thermal boundary-layer thickness is used. 

At the stagnation point of any blunt-nosed cylindrical body, con- 
ditions are the same as those at the stagnation point of a plate normal 
to the flow. Therefore m m 1, but the value of 6* is unknown. 

However, there exists at the stagnation point a unique relation 
v^5* = F^(du s */dx*) J which may, for instance, be read along the 

abscissa in figure 3 or in figure 5. Squaring this equation and dividing 
both sides by (du*/cLx*)6* 2 result in 


These relations are presented in figure 7 for the dimensionless convec- 
tion boundary -layer thickness and in figure 8 for the dimensionless 
thermal boundary -layer thickness. 

By use of these charts, the calculations for any profile can be made 
in a very simple manner for either the dimensionless convection or the 
dimensionless thermal boundary- layer thickness* The method of solution 
for the convection thicknesses is described subsequently. For the 
thermal thickness, the procedure is the same. 

The values of u s and du s /dx must be found for the cylinder pro- 
file under consideration either by measurement or by a solution of the 
inviscid-f low equations. The coolant velocity is prescribed by the 
porosity of the wall and by the pressure distribution around the profile. 
From these terms, the values of u*, du*/dx*, and v* can be calcu- 

culated. The value of 5* at the stagnation point can be determined 



dx* dx * 
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v*2 

from figure 7 in the following way: The value of — — - — - is computed, 

dug/dx 

and the corresponding value of (du*/dx*) is read from figure 7. A 

simple algebraic operation then yields the desired value of 6 * at the 
stagnation point. 

)|( 

The dimensionless convection boundary-layer thickness 6 C along 
the cylindrical surface is determined from equation (33 ) } for the 
numerical evaluation presented herein, this equation was solved by the 
method of isoclines with the aid of figure 3, depending upon which ratio 
of stream to wall temperature applied. Equation (33) determines the 
direction of the tangents to the different 5 C -curves which satisfy the 
equation. The task is to find that curve which contains the 5 c -value ^ 
previously calculated for the stagnation point. For chosen values of x 
and 5*, values of (du*/dx*) 5 C * 2 and are computed and the 

value of M is read from the appropriate part of figure 3. Equation (33) 
then gives the slope of the tangent at this selected value of ^ x for 
the assumed &*. Several values of 6* are used for this x . The 

same calculations are repeated for other values of x . If the distance 
between these x*- values is chosen small enough, an accurate curve of 
6* against x* can be drawn which starts at the desired previously 
calculated value of 5* at the stagnation point and which will have the 
correct slope at each value of x* considered. Figure 9 illustrates 
this method of solution. Values of N can then be obtained for each of 
the correct &q— values and the considered v-^,- value for each x from 
figure 4 after (du*/dx*) 5 C * 2 and v^S* are computed (the ratio of 

stream to wall temperature under consideration determines which part of 
figure 4 should be used) . The value of Nu/ a/Rcq can finally be 
obtained from equation (35). 

The same calculation procedure can be used when the dimensionless 
thermal boundary-layer thickness is considered. Figure 8 is used for 
the determination of the value of 5* at the stagnation point j figure 5 
is used to determine and figure 6 is used to determine N* The par- 

'kpQuqoj'* ratio of stream to wall temperature under consideration deter- 
mines which parts of these figures apply for the calculation of the 
values of M and N* Finally, equation (35) gives the desired value of 

Nu/ VRe 0 - 


Charts and Calculation Procedure for Prescribed Wall Temperature 

The heat-transfer coefficients determined by the values of Nu/ a/f^q 
can now be used to calculate the surface temperature of the cylinder 
when the outside stream temperature and the temperature with which the 
coolant is supplied to the interior of the cylinder are known. For this 
purpose a heat balance for an element of the wall as shown in figure 10 
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is set up. The cylindrical volume element considered may have two plane 
surfaces, one surface (l) coinciding with the outside surface of the 
cylinder wall and the other (2) apart from the inside surface of the wall 
hy such a distance that it is situated outside the boundary layer present 
on this side. (The inside surface has to be considered as a surface of 
a wall to which suction is applied and on which a boundary layer builds 
up as shown in reference 29.) The mantle surface (3) of the cylinder may 
be normal to the wall surfaces. Heat is carried by convection with the 
cooling air through surfaces 1 and 2. The amount per unit time is indi- 
cated in figure 10. It is assumed that the coolant is heated up to the 
wall surface temperature T w when it leaves the wall. This assumption 
is usually well fulfilled. Heat .will he also transferred hy conduction 
through the fluid layers immediately adjacent to the outside wall sur- 


faces, the amount being 




dA. In addition, heat may he trans- 


ferred to the outside wall hy radiation; it may he q r dA. Heat may 
also flow into the volume element hy conduction in the solid material or 
hy transverse flow of the cooling air. The sum of all these individual 
flows may he dA. Then the heat balance is 


q. r + 0^. + c pPw v w^c “ c pPw v w^w “ ^w^By 


w 


The heat -k w (dT/dy) w transferred per unit area from the gas to the 
wall was in this report expressed hy a heat-transfer coefficient 


h(^-T S > - 


Combining these two equations results in 

*r + <3c + h ( T s- T w) = c pPw v w( T w“ T c) 


(38) 


This equation permits a calculation of the wall temperature for any 
place on the cylindrical surface when the coolant velocity v w is pre- 
scribed, when the local radiative heat transfer q r and the conductive 
heat flow q c are known, and when the heat -transfer coefficient h has 
been obtained. The conductive heat flow q c is usually small and can 
be neglected. Such a calculation results in a wall surface temperature 
which generally will vary along the circumference of the cylinder. When 
the variations are large, the temperature distribution obtained can be 
regarded only as an approximation, since the wedge solutions (refer- 
ences 22 and 25) on which the method in this paper is based were obtained 
for the case of a constant wall temperature. 

Usually, however, the problem which faces the designer in an appli- 
cation is somewhat different from the one treated. The purpose of 
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"transpiration cooling is mostly to keep the wall temperature of some 
structural element below the limits which the material can withstand.. 

On the other hand., the amount of coolant almost always must be kept 
small, which means that local overcooling should be avoided. For the 
wall surface under these conditions, a temperature is prescribed which 
should be uniform about the circumference of the cylinder and the prob- 
lem is to find that distribution of the coolant velocity v w which 
results in the desired wall temperature. Generally, such an investiga- 
tion requires a trial-and-error procedure which is very involved. The 
procedure becomes simple and straightforward, however, when the radiative 
heat flow q r and the conductive heat flow q c can be neglected. Such 
a solution then is useful also as a starting point for the trial-and- 
error procedure for the case when radiation is present. 

The heat balance (equation (38)) can be transformed to 

. p v T -T 

h _ H w w w c 

c pPs u s,0 Ps u s,0 T s“ T w 

when q = q c = 0- The ratio of temperature differences in this equa- 
tion is r now a prescribed value. A similar ratio (T S -T W )/(T S -T C ) often 

appears in turbine-cooling work and is denoted by cp. Introduction of 
this value and conversion to dimensionless values results in 

-JH- = Pr ilSP (40) 

W 'P 

Another expression for Nu/ n/Rcq is given by equation (35) • Combining 
both equations gives 

N = v^6* Pr ~ (41) 

j|[ 

This equation expresses a relation between the parameters N and v^8 
in figures 4 and 6 which may be used to insert lines of constant cp 
into these figures* With the use of these lines the calculation pro- 
cedure for any specific problem becomes quite simple. The procedure 
will be described for T s /T w =1 (or near 1) and with the use of the 

convection boundary-layer thickness 5 C - The prescribed temperatures 
fix the value of cp. 

At the stagnation point, m = 1 and du*/dx* is known. In fig- 
ure 4(a) the intersection between the line m = 1 and the line for the 
prescribed cp determines v^S * and (du*/dx*) 6* and from both 

values, 5* and v* may be calculated. 
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The method of isoclines may again he used to determine the develop- 
ment of the boundary layer along the cylindrical surface. The use of 
this method implies that the gradient db*/dx* has to be determined for 

any pair of values x* and &£. For an assumed 5 *, the value v^&£ 
can be found in figure 4(a) as the value on the prescribed cp -curve above 
the known abscissa value (du^/dx*) 5^. Figure 3(a) then gives M 
and equation (33), the gradient db^/dx*. A plot similar to figure 9 
determines the boundary -layer thickness, and the values v belonging 
to these boundary-layer thicknesses represent the coolant-flow distri- 
bution for the particular temperature- difference ratio cp . 


NUMERICAL EVALUATIONS AND COMPARISONS WITH KNOWN RESULTS 

Solid Surfaces 

The results of the outlined procedure to calculate local heat- 
transfer coefficients have to be compared with experimental results or 
calculations by some other method in order to check the accuracy. The 
only cylindrical shape for which experimental data or solutions of the 
boundary -layer equations suitable for such a comparison are available 
seems to be the cylinder with a circular cross section. Accordingly, 
local heat -transfer coefficients were calculated by the method proposed 
in this report with the use of the dimenionless thermal boundary- layer 
thickness as well as of the dimensionless convection boundary- layer 
thickness. The results of these calculations are plotted in figure 11 
over the dimensionless distance from the stagnation point. Also 
inserted in the figure is a curve representing the average curve through 
the experimentally determined local heat-transfer coefficients mentioned 
in reference 30. It was shown in reference 22 that the measurements 
correlated well into a single curve when the experiments with very high 
Reynolds numbers near the critical value for transition to turbulence 
within the boundary layer were excluded. The tests with high Reynolds 
numbers gave values of Nu/ a/R^q over the whole upstream side 

of the cylinder were about 10 percent higher than the ones for the lower 
Reynolds numbers. The same behavior was found in references 31 and 32 
in which it is shown that an increase up to 50 percent in the heat- 
transfer coefficients over the expected laminar values was caused by the 
turbulence level in the wind tunnels used. The result of a solution of 
the boundary-layer equation as presented in reference 1 is also included 
in figure 11. This method solves the boundary-layer equations and 
obtains results as a series in the distance along the surface. Also 
inserted are values obtained by use of the Pohlhausen flat-plate solu- 
tion when the free-stream velocity is based on the local values and 
results obtained by the methods of references 12 and 13. Heat-transfer 
coefficients on wedges with the same local stream velocity and velocity 
gradient at the same distance from the stagnation point are also included 
(reference 8) . Appendix C explains how these wedge solutions were 
obtained. 
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On a cylinder with a circular cross section, separation occurs in 
the subcritical range near the value x* = 0.7. The stream velocity 
distribution around the surface of the cylinder which was needed for the 
calculations was obtained from pressure distributions given in refer- 
ence 30 and is contained in reference 22. 

It may be seen from figure 11 that the use of flat-plate values 
results in heat-transfer coefficients which are considerably lower than 
experimental values, whereas the methods in references 12 and 13 result 
in values which are too high. Much better agreement is found between 
the wedge heat-transfer coefficients and the experimental results, 
especially near the stagnation point. Farther downstream, the accuracy 
is improved by the method of this report. For the largest distance from 
the stagnation point, the use of the dimensionless thermal boundary- 
layer thickness results in values which are higher and the use of the 
dimensionless convection thickness, in values which are lower than the 
experimental ones. The values calculated by Frossling's solution of the 
boundary -layer equations are also higher than the experimental ones. 
Frossling's method has to be considered as an exact solution of the 
boundary- layer equations. In reference 22 it is recommended, on the 
basis of the good agreement between Frossling's curve and the values 
obtained by the use of the thermal boundary- layer thickness, that the 
method of the equivalent wedge flow be based on the thermal boundary- 
layer thickness. The values of the heat-transfer coefficients depend 
primarily on the velocity distribution in the stream around the cross 
section of the cylinder. The velocity distribution used for the calcu- 
lation on the circular cylinder is also shown on figure 11. The calcu- 
lations are made for a Prandtl number of 0.7, for a solid surface 
(v w = 0) and a temperature ratio T g /T w of 1, equivalent to the assump- 
tion of constant property values. These calculations agree within 
5 percent with the exact calculation and within 8 percent with experi- 
ment when the immediate neighborhood of the separation point is excluded. 
Similar comparisons have already been made in reference 2 for a gas with 
a Prandtl number of 1 and a different velocity distribution (see fig. 1). 
This comparison shows that the method proposed by Squire (reference 16) 
gives heat -transfer coefficients which agree with the exact boundary- 
layer solution to about the same degree as those of the method of the 
equivalent wedge flow. The same fact holds for the method indicated in 
references 15 and 17 especially with the improvement given in refer- 
ence 4. It can be stated in summary, therefore, that a number of methods 
exist today which at least for the circular cylinder permit the deter- 
mination of heat-transfer coefficients on solid surfaces in the laminar 
region of a gas having constant property values with a very good accuracy. 
The advantage of the equivalent wedge flow method over those methods 
just discussed is that it gives solutions in a very short time and that 
it can be readily extended to include variable property values and 
transpiration cooling, as was done in this report. The wedge solution, 
according to reference 7, is still more rapid} however, the results 
differ from the experimental values up to 15 percent. 
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Figure 12 gives the analogous results for an elliptic cylinder with 
the axis ratio 1:2. It may be observed that heat-transfer coefficients 
on wedges differ only slightly from those obtained for equivalent wedge - 
type flow, whereas the flat-plate values and the ones calculated with 
references 12 and 13 are considerably different. No experimental 
results nor solutions of the boundary -layer equations for a cylinder 
with such a cross section which could be compared with the approximate 
solutions are known to the authors. A calculation with Krou.jiline ' s 
method presented in reference 22 agrees well with the solutions obtained 
with the equivalent wedge -type flow method. Separation of the flow 
occurs on such a profile near x* = 0.8. The stream velocities used are 
calculated values contained in reference 22. 

The agreement between the wedge solutions and the results obtained 
by the method herein is still closer for the elliptic cylinder with axis 
ratio 1:4 as can be seen from figure 13- The reason for this fact is 
the type of stream velocity variation occurring on elliptic cylinders. 
Flow separation occurs on this cylinder near x* = 0.85. The curves in 
figures 12 and 13 show that the stream velocity is comparatively con- 
stant over a considerable part of its circumference after a steep 
increase near the stagnation point. This behavior is the more pro- 
nounced for an axis ratio of 1:4 than for one of 1:2. An inspection of 
figure 13 shows that, apart from the region near the stagnation point, 
even the flat-plate values give a reasonably good approximation. Cal- 
culations obtained by use of the dimensionless thermal boundary-layer 
thickness extended to the flow separation point, whereas those for the 
dimensionless convection boundary-layer thickness did not. It therefore 
appears advisable to use the dimensionless thermal boundary- layer 
thickness. 

Experimental heat-transfer coefficients found at the University of 
California for an elliptic cylinder with an axis ratio of 1:4 (ref- 
erence 33) are about 50 percent higher than the theoretical values 
shown in figure 13. There are several reasons for this discrepancy. 

The measured stream velocity distribution was different from the one on 
which the present calculations are based, probably because of a limited 
width of the wind tunnel. The cylinder in the experimental investiga- 
tion was heated by an electric resistance which produced a constant heat 
flow through the surface per unit area. Accordingly, the surface tem- 
perature varied along the circumference of the cylinder, being lowest 
at the forward stagnation point and increasing in downstream direction. 
Calculations in reference 33 indicate that the higher values found in 
the tests are mostly due to this fact. Another increase of the experi- 
mental heat -transfer coefficients may again be connected with the tur- 
bulence level in the wind tunnel used as discussed in connection with 
the test results on circular cylinders. 


22 


NACA TN 2733 


From figures 11 to 13, it may be concluded that, for cylinders with 
a stream velocity which is fairly constant over the greater part of the 
circumference, local heat-transfer coefficients may be obtained with 
good accuracy from wedge solutions. In the region in which the stream 
velocity variation is considerable, the method of the equivalent wedge 
flow gives heat-transfer coefficients with an accuracy sufficient for 
engineering purposes. 


Porous Surfaces 

Heat -transfer coefficients were calculated by the method of the 
equivalent wedge flow for cylinders with circular and elliptic cross 
sections for transpiration- cooled surfaces and different temperature 
ratios T s /T w by using either the thermal or the convection boundary- 

layer thickness (figs. 14 to 18). The use of both boundary- layer thick- 
nesses gives different results only for large distances from the stag- 
nation point. The variation of the heat-transfer coefficients with the 
ratio of stream to wall temperature is comparatively small for solid 
surfaces. This result is in agreement with previous findings. For 
transpiration-cooled surfaces, however, the effect of the temperature 
ratio on the heat -transfer coefficients becomes more pronounced, espe- 
cially on cylinders with nearly circular cross sections. In refer- 
ence 24 the case of transpiration cooling with small temperature differ- 
ences is calculated; this reference includes the effect of the tempera- 
ture ratio by a correction factor which is based on the assumption that 
this effect is the same as that determined experimentally for impermeable 
surfaces. A comparison of results shows that the procedure in refer- 
ence 24 underestimates the effect of temperature ratio for transpiration- 
cooled surfaces. In addition, it can be observed that transpiration 
cooling results in a considerable decrease of the heat-transfer coeffi- 
cients. A larger amount of coolant flow is necessary to reduce the 
heat-transfer coefficients by the same amounts in regions in which the 
heat-transfer coefficients are large. Such a region exists at the stag- 
nation point on the cylinder with the axis ratio 1:2, and especially on 
the cylinder with the axis ratio 1:4. 


The variation in coolant flow required to maintain constant wall 
temperature for transpiration-cooled cylinders with circular and 
elliptic cross sections is shown in figure 19. The calculations were 
made for a temperature ratio T s /T w of 1, a value of cp of 0.5, and 
a Prandtl number Pr of 0.7. Figure 19 shows that the highest local 
coolant -flow rates are necessary near the stagnation point in order to 
keep the wall temperature down at that place- The magnitude of the 
coolant-flow rate at the stagnation point is proportional to the square 
root of the velocity gradient du^dx ; this in turn is determined 
mainly by the value of the radius of curvature at this point. As this 
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radius of curvature decreases, the required coolant flow increases. 

This is in agreement with figure 19, which shows that the maximum coolant 
flow is required at the stagnation point of the elliptic cylinder with 
the 1:4 axis ratio. Downstream of the stagnation point, the flow rates 
decrease for each cylinder. Figure 19 also shows that the use of the 
thermal rather than the convection boundary -layer thickness results in 
only a very minor increase in coolant flow required to maintain the 
circular cylinder wall at a constant temperature. 


EXTENSION OF CALCULATION TO HIGH-VELOCITY FLOW 

The heat generated by internal friction was neglected in equa- 
tion (5) according to the assumption of small velocities. The equation 

q = h(T g -T w ) (42) 

gives the heat-transfer coefficient for this case. It was already 
explained that the inclusion of the internal friction for a gas with a 
Prandtl number of 1 results only in the change that the temperature T 
in equation (5) and the temperature T s in equation (42) are now total 
temperatures, as long as the property values my be regarded constant. 

The heat-transfer coefficients determined in this report my be used in 
this case. It was shown in reference 34 with the use of results obtained 
in reference 35 that the heat-transfer coefficients determined for low- 
velocity flow apply to high-velocity flow up to a Mach number of about 4 
for a gas with a Prandtl number different from 1, when the stream veloc- 
ity is constant (flat plate) and the heat flow is not too large. The 
heat-transfer coefficient, however, has now to be defined by the 
equation 


H = MT ad -T w ) (43) 

in which the temperature T ad denotes the value which an unheated plate 
assumes in the high-velocity flow. The adiabatic wall temperature my 
be determined from the recovery factor 


r 0 = 


T -,-T 
ad s 


c T,s" t s 


(44) 


which was found to be equal to /^/Pr for laminar flow and for Prandtl 
numbers not too far from 1. The difference between the total and the 
static temperatures in the stream is connected with the stream velocity 
by the equation 

2 


u c 


^T,s _ ^s " 


2c t 


(45) 
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For the flat plate with a constant stream velocity, the adiahatic wall 
temperature is therefore constant. 

Conditions are more involved on a cylinder with a stream velocity 
which varies along its circumference. Even when the recovery factor is 
assumed to he constant, equations (44) and (45) give an adiahatic wall 
temperature which varies along the circumference of the cylinder. The 
fact that the low-velocity heat-transfer coefficients also represented 
the high-velocity values on a flat plate, however, followed from the 
fact that the energy equation for constant property values is linear in 
T, and that a general solution of the nonhomogene ous equation describing 
the heat transfer including the internal friction could therefore he 
obtained by superposition of the solution of the homogeneous equation 
valid for small velocities and a particular solution of the nonhomo - 
geneous equation. Such a superposition results in a constant wall tem- 
perature on the flat plate when the solution of the homogeneous equation 
for constant wall temperature and the one describing the adiabatic wall 
temperature is used, since the adiabatic wall temperature is also con- 
stant. For a cylinder with an arbitrary cross section, however, the 
adiabatic wall temperature which represents a particular solution of the 
nonhomogene ous equation varies along the circumference. Therefore, a 
superposition of this particular solution with the low-velocity solu- 
tions for constant wall temperature does not give a constant wall tem- 
perature, which was specified for the problems investigated in this 
report. Accordingly, the heat transfer has now to be calculated with 
the equation 


q. = k(T eff -T w ) (46) 

in which T e ff has to be determined for constant wall temperature con- 
ditions, namely, as the temperature which a particular spot along the 
surface for which the heat-transfer coefficient is to be determined 
assumes when the heat flow through the wall at this particular spot is 
zero and the wall temperature along the circumference of the cylinder is 
constant. 

For flow around wedges, this temperature, which may be referred to 
as the T1 effective temperature,” can be found from the results in refer- 
ence 23. It is also determined for several cases in reference 36. The 
calculation procedure which determines this effective wall temperature 
from reference 23 is described in appendix D* The calculation shows 
that this temperature may be again expressed by a recovery factor 

T eff _T s 

r °°" t t , s -t s 


( 47 ) 
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The index oo is used to indicate that such a recovery factor could be 
determined experimentally by a model made of a material with a very 
large heat conductivity so that the internal heat conduction would elim- 
inate all temperature differences along the surface. On the other hand, 
the recovery factor describing the adiabatic wall temperature in equa- 
tion (44) has to be determined experimentally by a model made of a 
material with an infinitely small heat conductivity so as to eliminate 
internal heat flow. Values for the recovery factor r^ determining 
the effective temperature of a wedge are presented in figure 20. The 
recovery factors Tq describing the adiabatic wall temperature accord- 
ing to equation (44) have been calculated for wedges in reference 7. 

This calculation had resulted in values which decreased slightly with 
increas ing Euler number m. Repetition of these calculations on an 
electric computing machine, however, according to a communication from 
Arthur N. Tifford of Ohio State University, showed that the recovery 
factors for the adiabatic wall temperature are practically independent 
of the Euler number and have the same values as the recovery factor r^ 
shown in figure 20 for an Euler number m equal to zero. 

The consideration up to now dealt with solid surfaces. No infor- 
mation was found in the literature on recovery factors for transpiration- 
cooled surfaces. Some recovery factors were therefore determined for a 
transpiration-cooled flat plate and a flow with constant property values 
(the same for outside and coolant flow) by an integration of the 
boundary- layer equation (4) . The integration was carried out in the 
same way as in reference 5. The dimensionless stream function f and 
its second derivative were taken from reference 29. The results of 
this calculation are presented in figure 21 and the following table 
where T s /t w - 1 and Er = 0.7: 



Recovery factor 

-1 

0.713 

-.75 

.750 

-.50 

.786 

0 

.838 

.50 

.874 

1 

.900 


The figure shows that the recovery factors decrease considerably with 
increasing coolant flow. The calculations were extended to positive 
values of f w which apply to a surface with suction. 

It might be worthwhile to mention that the accurate determination 
of the adiabatic or effective wall temperature appreciably influences 
the heat flow as calculated by equation (43) only when the difference 
Tad-Tw is °f the same order of magnitude as or of a small order of 
magnitude than the difference T^ s -T g (see also appendix D) • 
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RESULTS AND CONCLUSIONS 

An approximate method for "the calculation of heat transfer in the 
laminar region around cylinders of arbitrary cross section was presented. 
The method, called the equivalent wedge -type flow method, is based on 
exact solutions of the laminar boundary-layer equations for wedge-type 
flow and takes into account the influence of large temperature differ- 
ences between the flow and the cylinder wall and the influence of tran- 
spiration cooling. The use of prepared charts reduces calculations to 
a graphical solution of an ordinary first-order differential equation. 

The method can be based either on the convection thickness or on the 
thermal thickness of the boundary layer. The results of calculations 
based on one thickness differ slightly from those based on the other 
thickness. There are not enough experimental data available to decide 
which boundary- layer thickness should be used. Near the separation 
point, however, the results obtained with the thermal boundary -layer 
thickness seem somewhat more plausible- 

The method was applied to circular and elliptic cylinders. The 
following results and conclusions are given: 

X. Results of experiments and exact calculations were available 
only for circular cylinders with solid surfaces. Calculations based on 
the present method and on the thermal boundary-layer thickness agreed 
within 5 percent with the exact calculation and within 8 percent with 
experiment when the immediate neighborhood of the separation point was 
excluded. 

2. With the present method, heat-transfer coefficients may be 
obtained without a knowledge of the flow boundary layer. Consequently, 
such calculations are more rapid than those based on the momentum and 
heat-flow equations. 

3. Heat-transfer coefficients determined from wedge solutions 
agreed on the circular cylinder within 15 percent with the results of 
experiments. The calculation procedure is still more rapid. 

4. For elliptic cylinders, the differences between the results of 
calculations with the various methods decreased as the axis ratio 
increased from 1:2 to 1:4. 

5. The development of the boundary layer is determined by the 
velocity distribution around the cylinder. The accuracy which has to be 
expected for the results of calculations with the different methods 
will therefore depend on the character of the velocity distribution. 

6- For cylinders with solid walls, the variation of the heat- 
transfer coefficients with ratio of stream to wall temperature was 
comparatively small. 
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7. For transpiration-cooled surfaces, the effect of temperature 
ratio on heat-transfer coefficients became pronounced, especially on 
cylinders with nearly circular cross sections. 

8. A considerable decrease in heat-transfer coefficients accompanied 
transpiration cooling. 

9. The influence of transpiration cooling on the recovery factor 
was investigated for a flat plate and constant property values. It was 
found that the recovery factor decreased considerably with increasing 
coolant flow. 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, March 19, 1952 
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APPENDIX A 
SYMBOLS 

The following symbols are used in this report: 
dimensionless wall temperature gradient taken from references 22 

“ d 23 > Mi % 

constant 

specific heat at constant pressure 
function 

dimensionless stream function, (p^tW h w P w xu g 

heat-transfer coefficient 
thermal conductivity 

characteristic dimension (major axis of cylinder) 

— i n b 2 = m ( 5* 2 , v*5*j (see equation (34)) 

Euler number, u s = ^ xm 

( du * \ 

= ^-^-*,6*2, v w & / ( see e( i uatiori ( 37 )) 

Nusselt number, hL/ky. 

Prandtl number, c^p/k 

pressure 

heat flow 

approximated heat flow 
heat flow by conduction 
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1 r 

Reo 

r O 

r co 

T 

T s 

u 

u s 

u s 

V 

V* 

X 

X* 

y 

z 


p 


5 

5 c 

5? 


heat flow by radiation 
Reynolds number, u s ,O l Po/mw 

recovery factor defined by (T a ^-T s )/(T iji^ s -T s ) (equation (44)) 

recovery factor defined by (T e ff -T s )/(Tt,s“T s ) (equation (47)) 

temperature in boundary layer 

temperature in stream 

velocity component along surface 

free-stream velocity 

dimensionless mass velocity in free stream, P w u s/ p 0 u s,0 
velocity component normal to surface 

P w v 

dimensionless velocity normal to surface, a P e n 

P0 u s,0 * ° 

distance from stagnation point along surface 

dimensionless distance from stagnation point along surface, x/L 
distance normal to surface 

dimensionless boundary-layer coordinate taken from references 22 
and 23, iJ~ 2 ~ 1 

1 IS 

/\[2lp x N v 

pressure gradient parameter, 2m/(m+l) 
boundary-layer thickness 

dimensionless boundary-layer thickness, (8 /l) /\j ReQ 
convection boundary-layer thickness (equation (15)) 
dimensionless convection boundary-layer thickness, (B c /l) AjReo 
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displacement boundary-layer thickness (equation (13)) 


8^ momentum boundary-layer thickness (equation (14)) 

8 t thermal boundary-layer thickness (equation (16)) 

8* dimensionless thermal boundary-layer thickness, (8 -^/L) a/Rcq 
q dimensionless boundary-layer coordinate, y ; 


% 

9 


Pw u s 


|Pw u S 




T-T. 


w 


dimensionless temperature -difference ratio, - — — 

i s _i v 


d dimensionless temperature -difference ratio. 


T-T, 


T T,s“ T s 


m+1 


A dimensionless stream function taken from reference 26 * 


w 


\i absolute viscosity 

V kinematic viscosity, p/ p 

£ distance along wedge, taken from references 22 and 23 
p density 

T -T 
± s w 

cp dimensionless temperature -difference ratio, 

1 s” i c 

^ stream function 

Subscripts : 
ad adiabatic 

c coolant, when used with T 

eff effective 
s stream * 


[ 
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T total 

w wall 

0 except when used with r, refers to a fixed point in the stream 
Superscripts : 

m exponent of distance along surface from stagnation point for stream 
velocity, u s = Cx 111 

' denotes differentiation with respect to t] 
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APPENDIX B 

EVALUATION OF HEAT -FLOW EQUATION 

The energy equation (5) will be integrated along y throughout the 
boundary layer under the conditions of small Mach number , constant wall 
temperature, and constant specific heat 

c pJ 0 ( pu s + 0V f) dy= J^ ^( k f) dy 

The first term on the left side can be transformed by partial differen- 
tiation to 


pu g (PUT) - T 


An analogous transformation of the second term and consideration that 
the temperature gradient dT/dy is zero outside the boundary layer 
(for y = co ) result in 


w 


r® 

P 00 

00 

^ (puT)dy - | 

o x J 

T d(pu) dy + pvT 
ox 

0 

0 o 


0 


T 5 (P V ) dy 

oy 


Mar) 

C p W / w 


The second and fourth terms cancel because of the continuity equa- 
tion (3). In the first term, the sequence of differentiation and inte- 
gration can be reversed. Introduction of the convection thickness of 
the boundary layer leads finally to the integrated heat -flow equation. 


d , R v _ ^w (be \ 

p s dx ^ u s 5 c) ‘ p w y w C p \3y/ w 


(Bl) 


It will now be proved that equation (26), used for the method of the 
equivalent wedge-type flow, is the same as this integrated heat-flow 
equation when the convection thickness for the boundary layer is used. 
Equation (Bl) may be transformed by partial differentiation of the 
first term into 


d&c 

p s u s dF 


du c 

+ ^ ar 


- Pw v w 


= ^w ( as | 
c T) W / 


(B2) 
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For wedge -type flow, the convection thickness is given by the expression 


Introducing this expression as well as equations (9) and (12) into 
equation (B2) gives the equation 


which interconnects the convection thickness with the dimensionless 
temperature gradient at the wall. The gradient of the convection thick- 
ness may now be determined from the integrated energy equation (B2) when 
the expressions in this equation are transformed to the new variables 


Replacing the nondimens ional temperature gradient in this equation by 
equation (B5) results in 


which is the same as equation (26). 

It can also be proved by a completely analogous calculation that 
the method of the equivalent wedge -type flow, when it is used to cal- 
culate the momentum thickness of the flow boundary layer, satisfies the 
integrated momentum equation which is obtained from equation (2) by an 
integration over y in a manner similar to the derivation of equa- 



(B3) 


Differentiation of this equation gives 



(B4) 



(B5) 


C A 

- — = -m T) _ A 

dx *c 




tion (Bl) 



w 


34 


NACA TN 2733 


APPENDIX C 

DETERMINATION OF WEDGE SOLUTIONS 

The wedge solutions which were used as a first approximation in 
figures 11 to 13 can he obtained very easily with the use of figure 22 
reproduced from reference 9. The heat-transfer coefficient has to be 
determined on a wedge which has the same stream velocity and its gradient 
at the same distance from the stagnation point as the real profile. The 
Euler number for this wedge can be found from equation (23). In the 
dimensionless coordinates it is 

du 

y S 

dx 


The parameter f w which determines the coolant flow through the porous 

wall is found from equation (20), which reads, when converted to dimen- 
sionless quantities, 


* da 
s dx 


(Cl) 



(C2) 


The value 
Nu/^Re o is 


(Nu///Re^) ^ can he determined 

finally obtained by multiplication 


from figure 22 , 


and 


When the temperature ratio cp is prescribed, figure 23 reproduced 
from reference 9 can be used to obtain the parameter f w for any Euler 
number m. Equation (C2) then determines the value v* and the dis- 
tribution of the required coolant flow along the profile. 
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APPENDIX D 


DETERMINATION OF EFFECTIVE WALL TEMPERATURE 


It is shown in reference 23 that for high-velocity flow of a fluid 
with constant property values around a wedge with constant wall temper- 
ature, the temperature field can be expressed by the equation 


in which 8 represents the nondimensional temperature field for low- 
velocity flow and $, the nondimensional temperature field for high- 
velocity flow and a wall temperature equal to the total stream temper- 
ature. The heat flow from the wall, obtained by differentiating equa- 
tion (Dl), is 


With the transformations used in reference 23 (see also appendix E) 


T ~ ( T w" T T,s)( 1 ~ 0 ) + ( T T,s -T s) $ + T s 


(Dl) 



z 




y 



u = Cx 
s 


m 


(D3) 



equation (D2) can be transformed into 



This equation is to be brought into the form 



(D5) 
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A comparison of equations (D4) and (D5) gives 






(D6) 



from which the difference between the total and effective gas tempera- 
tures can be found. The expression 


defines this temperature difference and the recovery factor for the 
effective wall temperature. The nondimens ional temperature gradients 
appearing on the right side of this equation are presented in refer- 
ences 23 and 36. In this way, the values in figure 20 have been deter- 
mined. 

To obtain an estimate of the conditions under which the difference 
between the adiabatic wall temperature and the effective wall tempera- 
ture may be neglected, the heat flow into the wall will be approximated 
by the equation 


and the error of such an approximation will be determined. The ratio 
of the exact heat -flow equation (D2) to the one approximated by equa- 
tion (D8) is 


Introducing the recovery factor for the adiabatic wall temperature 




(D7) 




(D8) 



T T,s- T ad = (! - r oK T T,s- T s) 
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gives 


q.w _ 

~~ 


T T , s~ T s 

Tv~T a d 




CO 


For an Euler number equal to 1, which characterizes flow near a stagna- 
tion point and which, according to figure 20, shows a large difference 
between the recovery factors tq and r^, the error is smaller than 
5 percent when 


T 

T 


ad~ T w 

T,s'" T 


s 


is larger than 2.5. 
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APPENDIX E 

COMPARISON OF VARIABLES 

This appendix gives a comparison of the variables used in refer- 
ences 9 and 25 with the ones used in references 22, 23, and 26. All of 
these references deal with wedge -type flow 


Symbols from 

Symbols from 

Relation 

references 9 

references 22, 

among symbols 

and 25 

23, and 26 



m (Eu) 

p 

P = 

2m 

m+1 

f w 

X 

X = 

m+1 f 

- — x w 




. |m+l 

n 

z 

z = 


0* 

A 

A = 

J 2 0' 

w 



\m+l w 


The values used in this report are related to the ones in the pre- 
ceding references by the following equations: 
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Figure 2. - Sketch of cylinder indicating notation used. 
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(b) Tg/V 2. 

Figure 3. - Continued. Chart for use in determination of M for dimensionless convection boundary-layer thickness. Pr, 0.7. 
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(c) T s /T w , 4. 

Figure 3. - Concluded. Chart for use in determination of M for dimensionless convection boundary-layer thickness. Pr, 0.7. 
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Figure 4 . 


(a) Ta/T w , 1. 


Chart for use in determination of N 


for dimensionless convection boundary -layer thickness. 


Pr, 0.7. 
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(a) T s /T w , 1. 

Figure 5. - Chart for use in determination of M for dimensionless thermal boundary- layer thickness. Pr, 0.7. 
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Figure 5. - Continued. 
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Chart for use In determination of M for dimensionless thermal boundary- layer thickness. 
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(c) Tg/T w , 4. 

Figure 5. - Concluded. Chart for use In determination of M for dimensionless thermal boundary- layer thickness. Pr, 0.7. 
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Figure 6. 


(a) T a /T w , 1. 


Chart for use in determination of N 


for dimensionless thermal boundary- layer thickness. 


Pr, 0.7. 
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(o) Tg/T w , 4. 

Figure 6. - Concluded. Chart for use in determination of N for dimensionless thermal boundary- layer thickness. Pr, 0.7. 
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du. 


dx 


s R *2 
— 5 c 



Figure 7. - Chart for use in determination of dimensionless 
convection boundary-layer thickness at stagnation point. 


2447 



2447 


8S NACA TN 2733 


57 



* / * 

d.u s /dx 


Figure 8. - Chart for use in determination of dimensionless 
thermal boundary- layer thickness at stagnation point. 
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Figure 9. - Isocline solution of boundary -layer 

equation. 



Figure 10. - Cross section through part of cylinder wall 
used in setting up heat balance. 
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4 0 .1 *2 .3 * -4 .5 .6 .7 


x 

Figure 11. - Comparison of present method -with formerly used methods for calculation 
of local heat-transfer coefficients around circular cylinder with impermeable wall. 
P** y 0-7 y T g /T w , !• 
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Figure 12. - Comparison of present method with formerly used methods for calculation 
of local heat-transfer coefficients around elliptic cylinder with axis ratio of 
1:2 . 0$ JPr, 0 • 7 j T s /T w ^ 1* 
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Figure 13* - Comparison of present method with formerly used methods for calculation 
of local heat-transfer coefficients around elliptic cylinder with axis ratio of 
1:4. v w , 0; Pr, 0.7; T s /T w , !• 
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Figure 14. - Local heat-transfer coefficients around transpiration- cooled circular 
cylinder determined hy use of dimensionless convection boundary -layer thickness for 
several cooling-air-flow rates and ratios of stream to wall temperature. 
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Figure 15. - Local heat-transfer coefficients around transpiration-cooled circular 
cylinder determined by use of dimensionless thermal boundary-layer thickness for 
several cooling-air-flow rates and ratios of stream to wall temperature. 
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Figure 16. - Local heat-transfer coefficients around transpiration-cooled elliptic 
cylinder with axis ratio of 1:2 determined by use of dimensionless convection 
boundary- layer thickness for several cooling-air-flow rates and ratios of stream 
to wall temperature. 
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Figure 17. - Local heat-transfer coefficients around transpiration-cooled elliptic 
cylinder with axis ratio of 1:2 determined by use of dimensionless thermal boundary- 
layer thickness for several cooling-air-flow rates and ratio of stream to wall 
temperature of 2. 
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Figure 18. - Local heat -transfer coefficients around transpiration- cooled elliptic 
cylinder with axis ratio of 1:4 determined hy use of dimensionless convection 
boundary- layer thickness for several cooling-air-flow rates and ratios of stream to 
wall temperature. 
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Figure 19. - Coolant flow required to maintain constant cylinder wall temperature. 

CP, 0.5; Pr, 0.7; T g /T w , 1. 
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Figure 20. - Recovery factor values determining the effective 

temperature of a wedge. 
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Figure 21. - Recovery factors for transpiration-cooled flat plate and constant 
fluid properties. Tg/T^, 1; Pr, 0.7. 
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Figure 22. - Heat transfer through laminar boundary layer with 
and without porous flow (reference 9). 
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